Quantum wires in magnetic field: A comparative study of the Hartree-Fock and the 

spin density functional approaches 
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We present a detailed comparison of the self-consistent calculations based on the Hartree-Fock 
and the spin density functional theory for a spit-gate quantum wire in the IQH regime. We demon- 
strate that both approaches provide qualitatively (and in most cases quantitatively) similar results 
for the spin-resolved electron density, spin polarization, spatial spin separation at the edges and the 
effective g factor. The both approach give the same values of the magnetic fields corresponding to 
the successive subband depopulation and qualitatively similar evolution of the magnetosubbands. 
Quantitatively, however, the HF and the DFT subbands are different (even though the correspond- 
ing total electron densities are practically the same). In contrast to the HF approach, the DFT 
calculations predict much larger spatial spin separation near the wire edge for the low magnetic 
fields (when the compressible strips for spinless electrons are not formed yet). In the opposite limit 
of the large fields, the Hatree-Fock and the DFT approaches give very similar values for the spatial 
spin separation. 

PACS numbers: 73.21. Hb, 73.43.Cd, 73.23.Ad 



Introduction. A detailed knowledge of energetics, spin 
splitting, magnetosubband and edge state structure in 
quantum wires is necessary for understanding and inter- 
pretation of a variety of magnetotransport phenomena 
in the integer quantum Hall (IQH) regime. A powerful 
tool to study electron-electron interaction and spin ef- 
fects in quantum wires is the mean-field approaches such 
as the Hartrec-Fock (HF) and the spin density functional 
theory (DFT)i. A number of studies based on these ap- 
proaches addressing various aspects of interacting elec- 
trons in the IQH regime in quantum wires have been 
reported recentlj^^^^ii^ia^^i^^ii^. However, in 
some cases different studies arrive to conflicting results 
and findings reported in some studies are not recovered 
in others. It is not clear whether the reason for such 
discrepancies is due to utilization of different approaches 
treating the exchange and correlation effects in different 
ways (i.e. HF vs spin DFT), or this difference is related 
to various approximations of different models (such as 
e.g. neglecting a global electrostatics, simplified models 
for screening, non self-consistent calculations, fixed filling 
factors, etc.). 

The aim of this Brief report is to resolve this issue by 
presenting a detailed comparison of the self-consistent 
calculations based on the HF method and the spin DFT 
approximation for a spit-gate quantum wire in the IQH 
regime. This includes a comparison of the magnetosub- 
band structure, electron densities, spin polarization and 
spatial spin separation as well as calculation of the effec- 
tive g factor. We stress that in our calculations we do 
not use any simplified assumptions concerning screening 
such that a global electrostatics of the system at hand 
is treated in an exact way. Note that comparative stud- 
ies of different approaches are common in treatment of 
electronic properties of quantum dots as they provide an 
important insight into the validity of applied methods 



and used approximations^^. At the same time, we are 
not aware of corresponding studies for the quantum wires 
in the IQH regime. Another motivation for the present 
study is that the DFT based approaches are often used 
for transport calculations. It has been argued recently 
that the standard DFT approaches might not be always 
suitable for this purpose (notably in the weak coupling 
regime), because of the inherent problems of the deriva- 
tive discontinuity problem and related uncompensated 
self-interaction errors in the DFT—. At the same time, 
the HF approach does not suffer from the above prob- 
lems and thus the utilization of the HF method in the 
transport calculations overcomes the limitations of the 
standard DFT. Thus, the detailed comparison of these 
two methods can provide a justification for the utiliza- 
tion of the HF instead of the DFT. 

Basics We consider an infinitely long split-gate 
GaAs/AlGaAs quantum wire in a perpendicular mag- 
netic field B where electrons are situated at the distance 
b below the surface. The HF equation for a single-particle 
wave function of spin cr, $^(r), readsi 



[Ho{r) + V,onf{y) + VH{y) + Vz] $^(r) 



(1) 



where r = {x,y), Ho{r) - 
is the kinetic energy in the 
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Landau gauge, 
TO* = 0.067me being the GaAs effective mass; a ~ 
describes spin-up and spin-down states, t i i- In the 
split-gate geometry the bare confining potential Vconf (y) 
due to the gates, donor layers and the Schottky barrier 
is well approximated by the parabolic confinement, 
Vconf (y) = Vo + ^ (woy)^, where Vq defines the bottom 
of the potential (we set the Fermi energy Ep = 0)^. 
The Zeeman energy is Vz = gj-itBa where ^ eh/2me 
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is the Bohr magneton and the bulk g factor of GaAs 
IS g = —0.44. The Hartree potential due to the elec- 
tron density n{y) = Tlia'^^^iv) (including the mirror 

charges) ^ Vniv) = -j^J dy'n{y')\n J^ff,,. ■ 
The non-local Fock operator is VFocfc(r,r') = 

~T^^I^\T.pf^^n^^)n*^^')^ ^J^C'^e the sum- 
mation is performed over all states and is the 

Fermi-Dirac distribution function. 

We assume the Blocli form of the wave function, 



(2) 



where k is the wave vector, ^n.k{yY describes the n-th 
transverse subband for the spin a. Substituting the Bloch 
function ^ into the HF equation ([T|) and integrating 
over the longitudinal coordinate x we arrive to the set of 
eigencquations for ip"^ ^ (y) ^ 



hk\ 



Vz 



/ ^Foc;c(2/,y')<fc(2/')rfy' = K,fc<,fc(2/), 

(3) 




where lDc is the cyclotron frequency and 
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FIG. 1: (Color online). (a),(b) Spatially resolved spin polar- 
ization of the electron density v) — . (c),(d) The number 
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with Kq being the modified Besscl function. Discrediting 
Eq. ([3]) we reduce the system of the integro-differential 
equations to the system of linear equations, which we 
solve numerically by standard methods in an iterative 
way until the self-consistent solution is achieved. Knowl- 
edge of the wave vectors /c^ for different states {n, cr} 
allows us to recover the subband structurei^, i.e. to cal- 
culate an overage position of the wave functions for 
different modes n, = hk!^JeB. 

Within the framework of the spin density functional 
theory, the Kohn-Sham equations for the single particle 
wave function ^'^(r) read 



[i/o(r) + y^(y)]<I>^(r)=£;^$-(r), 
V^iy) = V^onfiy) + Vniy) + Vz + 



(5) 
(6) 



where the first three terms in the effective confinement 
potential V'^{y) are the same as in the HF equation ((T|), 
and the last term corresponds to the exchange and cor- 
relation potential in the local spin density approxima- 
tion. It is given by the functional derivative V^^'' = 
{ne^c(n)} , where e^.^ (tt-) is the exchange and cor- 
relation energy functional and C(y) = "I T"! is the lo- 
cal spin polarization. 
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arrows indicate the magnetic fields corresponding to the mag- 
netosubband structure shown in Fig. (2] (e),(f) The effective 
g factor. ( g),(h) The spatial spin separation at the wire edge 
dsep (a definition of dsep is outlined in Fig. [2] (a) ). The 
left and right panels correspond respectively to the HF and 
the spin DFT approximations. The bare confining potential 
Vconfiy) is well approximated by the parabolic confinement 
with Vc) = — 85meV, hujc = 2meV, the distance to the surface 
b = 60nm, temperature T = IK. 



the wave functions the equations (O are solved sclf- 
consistently using the Green's function technique as de- 
scribed in detail in Rcf. [H (sec also Rcfs. 14]1T5|). 

Note that we find the self-consistent solutions for the 
DFT and the HF approaches using completely unrelated 
numerical methods. As a validity check we control that 
these different methods give identical results when we set 
^Fock ^^^^ ^xc' in respectively Eqs. ([3]) and © to zero 
and thus reduce both approaches to the standard spinless 
Hartree approximation (the latter is shown to reproduce 
welliii^ the Chklovskii et al)^ electrostatic treatment). 
Results and discussion. Figure [T] show a spatially re- 



solved difference in the electron density — and the 



total spin polarization P 
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calculated using the 
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AU the results presented below 
correspond to the parameterization of e^^ 
Tanatar and Ceperlej»i^. 



(n) given by 
Assuming the Bloch form of 
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HF and the spin DFT approaches for a representative 
quantum wire {n^jj = J n" {y)dy). A detailed analysis of 
the spin polarization in a split-gate quantum wire based 
on the spin-DFT approach is given in Refs. [islfl^ . For 
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the sake of comparison with the HF approximation we 
summarize below the main finding. The spin polarization 
is maximal for magnetic fields close to the depopulation 
of the even subbands, see Fig. [T](d). In this case the high- 
est occupied (odd) subband forms a compressible strip in 
the middle of the wire, and, therefore, the electron den- 
sity is mostly spin-polarized in the centre of the wire, 
see Fig. [5] (d). (We define the width of the compress- 
ible strip within the window \E ~ Ef\ ^ 27rfcT^iiiii^i^ 
which corresponds to the energy interval where the sub- 
bands are partially filled, < f''° < 1). Minima of 
the spin polarization correspond, instead, to depopula- 
tion of the odd subbands. At polarization minima the 
spin-up and spin-down subbands are fully (and practi- 
cally equally) occupied as their bottoms in the center of 
the wire arc situated below (or just on the border) of 
the energy window \E — Ep\ ^ 27rfcT, sec Fig. [2] (c). 
Because of this, the spin polarization in the centre of 
the wire is absent and it increases toward the edges of 
the wire because the spin-up and spin-down subbands 
intersect Ep at different distances from the wire center. 
For the case of spinlcss electrons the compressible strips 
are formed near the wire boundaries for sufficiently high 
magnetic field a^^i^^'^° . The exchange interaction, how- 
ever, completely or partially suppress the compressible 
strips leading to a spatial spin separation between the 
spin-up and spin-down states^^. This spatial separation 
causes a strong spin polarization near the boundaries 
which is clearly seen in Fig. [T] (b) for magnetic fields 
B > 2.75 T. This spin separation grows as the magnetic 
field increases because the width of the corresponding 
compressible strips for the spinless electons increases, see 
Figs. [11(b), (h). Note that the spin-DFT approach also 
predicts an almost constant (independent of B) spatial 
spin polarization even for lower fields B < 2.0 T, when 
the compressible strips are not formed yet. 

Let us now compare the results of the spin-DFT cal- 
culations with those based on the HF approach. It is 
remarkable that the both approaches give practically the 
same values of the magnetic fields corresponding to the 
successive subband depopulation, c.f. Figs. [l](c) and (d). 
Moreover, both approaches give practically the same to- 
tal electron density distribution n(y), see Fig. [2] (a) and 
(b). However, the spin-resolved densities are not always 
the same. In contrast to the spin-DFT approach, the HF 
calculation does not always exhibit a spatial spin polar- 
ization near the edges for the low fields (when the com- 
pressible strips for spinless electrons are not formed yet). 
This is the case for a quantum wire of Fig. [T]for B < 2.75 
T. The spin separation near the wire edges dsep is caused 
by the exchange interaction and it is known to depend 
on the steepness of the confinement potentia]~i2i: as the 
external confining potential becomes smoother, the spin 
separation growths. Figure [3] shows that while the spa- 
tial spin separation d^^'^ and d^Jp exhibit qualitatively 
same behavior as a function of the potential steepness, 
the DFT approaches predicts much larger spatial spin 
separation as compared to the HF method. Besides, the 
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FIG. 2: (Color online). (a),(b) The local filling factors v\v^ 
and V = = nh/eB) calculated within the HF 

and the spin DFT approaches for two representative magnetic 
fields (indicated by arrows in Fig. [l](c),(d). (c),(d) The mag- 
netosubband structure for spin-up and spin-down electrons 
calculated within the HF and the spin DFT approaches. The 
fat solid lines indicate the DFT effective confinement poten- 
tial , Eq. (|6} (note that because of the nonlocal character 
of the HF equations, it is not possible to define the effective 
confinement potential in the HF approach.) 



critical value of the potential steepness at which differ- 
ent spins become spatially separated is obviously lower 
in the HF approach. We stress that the difference be- 
tween d£p and d^Jp discussed above corresponds to the 
regime of the low fields, when the compressible strips for 
spinless electrons are not formed yet. For larger fields 
(corresponding to the formation of compressible strips 
for the spinlcss electrons), the Hatrce-Fock and the DFT 
approaches give very similar values for the spin separa- 
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FIG. 3: The spatial spin separation near the wire edge as 
a function of the confinement steepness calculated within the 
spin DFT and HF approaches in the regime of low fields [B = 
0.5T). The inset illustrates how the confinement steepness 
changes as the bottom of the parabolic confinement potential, 
Vb, varies {hujo is adjusted to keep the wire width constant, 
w = 500nm.) 
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tion, c.f. Figs. [l](g) and (h). In this case dsep is approx- 
imately equal to the width of the compressible strips for 
spinless electrons (see Ref. [3 for a detailed discussion 
of the suppression of the compressible strips by the ex- 
change interaction leading to the spatial spin polarization 
at the edges). Note that our study can not distinguish 
which approach gives a correct result for d^ep for the low 
field. This question can be resolved by a comparison to 
the exact results obtained by e.g. quantum Monte Carlo 
methods. We speculate at this point that it is the DFT 
approach that overestimates the spatial spin separation 
at lower fields. This conclusion is based on transport 
measurements on lateral quantum dots indicating that 
the spin polarized injection and detection by means of 
the spatial separation of spins can be achieved only in the 
edge state regime for sufficiently high magnetic field^^. 

Let us turn our attention to the subband structure. 
Figure [2] shows the subband structure for two representa- 
tive magnetic fields corresponding to the minimum and 
maximum of the spin polarization in a quantum wire. 
Qualitatively, the HF and the DFT subbands exhibit 
very similar features and evolve in a similar way as the 
magnetic field is changed. This includes the subband 
depopulation, the formation of the compressible strip in 
the middle of the wire and the subband separation at the 
edges. Quantitatively, however, the HF and the DFT 
subbands are different (even though the corresponding 
densities are practically the same, see Fig. [2](b), (d)). 
The most pronounced difference is that the consecutive 
subband separation for different spins in the DFT ap- 
proach is equal to fuvc, whereas the HF subband separa- 
tion exceeds this value. We attribute this difference to 
the nonlocal character of the exchange interaction in the 
HF approximation. Note that the HF subband separa- 
tion tends to hujc as the density increases because the 
exchange interaction becomes less pronounced for higher 
densities in comparison to the kinetic energy. 



Because the DFT and HF approaches give rather sim- 
ilar evolution of the magnetosubband structure, the cor- 
responding behavior of the total spin-polarization P and 
the effective g factor, g'^-^-^ , is also similar, see Figs. [1] 
(c), (d), and (e), (f). (We define the effective g fac- 
tor according to g"^^^ = (^{E\ f, — E^^ f^) / g^sB^ where 
the averaging is performed over the all k vectors and 
the occupied subbands n). The DFT approach gives 
a slightly higher value of P at the lower fields because 
of the enhanced spin polarization near the edges as dis- 
cussed above. The both approaches give quantitatively 
similar dependencies of g'^^^ as a function of magnetic 
field. Because g'^f^ is directly related to the subband 
spin splitting, the dependence oig'^f^ = g^^f{B) closely 
follows that of P = P{B), showing a well-known oscil- 
latory character with a periodicity of 1 /B related to the 
subband depopulation^^. A maximum value of g"^-^-^ w 15 
is reached close to magnetic fields corresponding to de- 
population of the even subbands, i.e. when the subband 
splitting in the wire center is maximal. 

Conclusion. We demonstrate that the spin DFT and 
the HF approaches provide qualitatively (and in most 
cases quantitatively) same description of a split-gate 
quantum wire in the IQH regime. This includes the elec- 
tron density, spin polarization and the effective g factor. 
The both approach give the same values of the magnetic 
fields corresponding to the successive subband depopula- 
tion and qualitatively similar evolution of the magneto- 
subbands. Quantitatively, however, the HF and the DFT 
subbands are different (even though the corresponding 
total electron densities are practically the same) . In con- 
trast to the HF approach, the DFT calculations predict 
much larger spatial spin separation near the wire edge for 
the low fields (when the compressible strips for spinless 
electrons are not formed yet). 
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